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Abstract 



Hadron production in lepton-nucleus interactions at high-energies is considered in 
framework of developing Monte Carlo event generator HARDPING (HARD Probe 
INteraction Generator). Such effects as formation length, energy loss and multiple 
^> . rescattering for produced hadrons and their constituents are implemented into the 

I HARDPING 2.0. Available data from HERMES collaboration on hadron production 

ly-^ ' in lepton-nucleus collisions are described by the present version of the HARDPING 

. generator in a reasonable agreement. 

O ■ Hadronisation of quarks and gluons is one of the most intriguing parts of nonperturbative 

^ ! QCD. Use of nuclear targets may allow to reveal important features of space-time picture of 

hadronisation, like hadron formation length and energy loss, see, e.g, for a review [1] [2] and 
references therein. The understanding of quark propagation in nuclear medium is crucial 
^ ■ for the interpretation of ultrarelativistic heavy ion collisions, as well as high energy proton- 

d ■ nucleus and lepton-nucleus interactions. To simplify interpretation of observable effects one 

can consider at the beginning hadron production in lepton scattering off nuclei. In case of 
deep inelastic scattering of lepton on nucleus there can be two stages of hadronisation. The 
first stage is predominantly perturbative. At this stage after hard scattering a struck a struck 
quark propagates through the nuclear medium being in point-like parton state experiencing a 
little attenuation only. This effect is known as Landau- Pomeranchuck-Migdal effect in QCD 
[31 m ini [6l [71 m [9l [TO] . At the end of the first stage, a pre-hadron state (a color dipole or 
constituent quark) is formed |Tll |13]. In the second stage pre-hadron state with smaller 
than hadron cross section interact with nuclear medium. There is finally formed hadron at 
the second stage. At the large enough energies of produced hadrons the nonperturbative 
stage of hadron formation is evolving beyond the nucleus 
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The aim of the work is to study these effects for case of lepton-nuclei colhsions using a 
developing Monte Carlo (MC) event generator. The generator HARDPING (HARD Probe 
INteraction Generator) is based on MC generators PYTHIA [H] and HIJING [15]. The 
first version of HARDPING describes experimental data on Drell-Yan reaction off nuclei 
reasonably well [Ml [IT]. It takes into account the effects related with interaction of pro- 
jectile hadron and its constituents in nuclear matter before hard scattering for lepton-pair 
production off nuclei. The second version of HARDPING, presented here, describes, in ad- 
dition, hadron production in lepton-nuclei interactions. It incorporates the following effects: 
formation length, energy loss and multiple soft rescatterings. 

The experimental results on semi-inclusive leptoproduction of hadrons off nuclei |18l [T9] 
are presented in terms of hadron multiplicity ratios R\j with nuclear [A) and deuteron [D) 
targets, as functions of virtual photon energy (i/), its fraction taken by hadron [zh] and 
hadron transverse momentum squared p]_: 



where N^^^ and N^^^ are yelds of inclusive deep-inelastic scattering leptons on nuclei A and 

yields of semi-inclusive hadrons as a function of x, here x is either Zh 
or In absence of nuclear effects, the ratio i?^ should be equal to 1. The experimental 
results show that this is the case at high transferred energy u [18] . 

It is well established from theoretical and experimental studies of hadron-nucleus colli- 
sions at high energy that hadrons are not produced at the point of collision but only after 
some "formation" length [Ij. In the Lund string fragmentation model, the production of 
hadrons is described as two stage process. At the first perturbative stage a pre-hadron at 
the end of the string is formed. On the next nonperturbative stage a hadron is formed. 
Before a pre-hadron is formed, the struck quark propagates through the nuclear matter with 
a very small cross section (in this work we neglect it). It takes some time at the pertur- 
bative stage to form a pre-hadron (formation time, tp or formation length, Ip). When the 
pre-hadron is formed, it interact with nuclear matter via pre-hadron cross section, which 
is different from hadronic cross section. And also it takes an extra time to form the final 
hadron from pre-hadron state. So, the formation length consist of two parts {Ip and /„), 
corresponding the two stages of hadronisation. 

There are two approaches to calculate formation length with two stages. The first is 
based on an oversimplified description of nonperturbative stage [T], while the second one [TT] 
is based on Lund string model, but neglecting the energy loss effects during the perturbative 
stage. The present work is based on the both above approaches with including the effect of 
energy loss at the perturbative stage and using Lund string model at the nonperturbative 
stage. 

In the first approach the distribution on pre-hadron formation time can be written in the 
following form [1]: 
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here tp is pre-hadron formation time, Zh is fraction of virtual photon energy carried out 
by the hadron, is virtual photon virtuality, v is virtual photon energy, Aqcd is QCD 
constant, i±) nieson wave function. Eg {tp) = u — AE (tp) is quark energy and AE (tp) 

is quark energy losses due to perturbation gluon radiation, 6 is delta-function, 5* {tp, z^jQ"^, v) 
is Sudakov suppression factor. AE {tp) can be taken in the following form: 



AE{t) = z/ / dej^^^}da 
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where 6 is step-function, as is strong coupling. This approach works only for leading hadrons 
with Zh > 0.5, (see FiglH the dashed line), which is not suitable for full MC simulation. 
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Figure 1: Formation length as a function of Zh and as a function of u. The dotted lines 
correspond to the first approach (B.Z. Kopeliovich et al. P) for Zh > 0.5, the dashed lines 
corresponds to the second approach (A. Accardi et al. [H]) and the solid lines corresponds 
to HARDPING 2.0 calculations 

The second approach is based on Lund string model for nonperturbative hadronisation 
neglecting its perturbative stage [H]. In this approach probability to have pre-hadron for- 
mation length Ip can be written in the following form: 

P {Ip] Zh, L) = i^z^ \— ^-^^ t1 X 



{5[lp-{l-zr,)L] + j^^e[{l-z,)L 

xe[ip]. 



(4) 
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where parameter C=0.3 [TT], and parameters k and L are string tension and ratio of the 
virtual photon energy to string tension L = u/k. 

In the presented here approach the effect of energy loss was incorporated into the HARD- 
PING 2.0 using PYTHIA MC implementation of parton shower for the perturbative stage 
and Lund string model for the nonperturbative one. On the FigH] the dependences of for- 
mation length on Zh and u are plotted, the dotted lines - the first approach [1], the dashed 
lines - the second approach [llj, the solid lines - HARDPING 2.0 simulation. 

During the perturbative stage, corresponding to the formation length Ip, a constituent 
quark (or pre-hadron) state is formed. It can interact with intranuclear nucleons via inelas- 
tic pre-hadronic cross-section (or inelastic quark- nucleon cross-section). At the end of the 
nonperturbative stage the observed hadron is formed. 

Produced pre-hadrons and hadrons can undergo soft collisions with intranuclear nucleons 
(with small momentum transfers: \t\ < 1 GeV'^). So, one has to take into account their soft 
multiple rescattering. 

The transverse momentum distribution of constituent quarks after one soft interaction 
can be parameterised in the following form [201 [ISl [IZ] : 

fp ipl) = ^e-^^- (5) 

where B = where (kp) is mean value of quark transverse momentum. fp{pi_) is a 

differential distribution of quark in quark-nucleon interaction normalised on unity. 

Probability to have no interactions between the points with coordinates (^z^ bj and (^z + X,b 
can be written in the next form: 

P{X-z,b)=e-<''^'') (6) 

where T (b, z, is: 

T{b,z,x) = {A-l)J^^^\{b,z')dz\ (7) 

p (b, z^ is nuclear density and a is quark-nucleon (pre-hadron-nucleon) or hadron-nucleon 
inelastic cross section. 

Simulations of lepton-nuclei collisions obtained by HARDPING 2.0 were compared with 
HERMES data [HI [19]. The results are shown on Figs. EH 

The performed simulations shown a reasonable agreement of MC model HARDPING 2.0 
with the experimental HERMES data [HI |T9]. This allowed to fix model parameters such 
as inelastic quark-nucleon (pre-hadron-nucleon) cross-section a = 10 mb and string tension 
k = 1.7 GeV/Fm. Comparision with EMC [21] and SLAC [221 123] data shown also a good 
agreement and it will be presented elsewhere. 

To summarise, the effects of the two-stage hadronisation and multiple soft interactions 
inside of nucleus for produced hadrons and their constituents were implemented into MC 
generator HARDPING 2.0. The developed MC generator HARDPING 2.0 is allowed to 
describe reasonably well the HERMES data [TH [19] on hadron production in positron-nucleus 
scattering at 27.6 GeV. 
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Figure 2: Multiplicity ratio {R^) of charged hadrons for krypton (Kr) and deuteron (D) 
targets as a function of and as a function of Zh at positron beam energy 27.6 GeV. 
The solid points correspond to HERMES data [18] and the open points are obtained by 
HARDPING 2.0 




Figure 3: Multiplicity ratio (-R^/) of vr"^ -mesons for xenon (Xe) and deuteron (D) targets 
as a function of p\ for different Zh values and as a function of Zh for different u values at 
positron beam energy 27.6 GeV. The solid points correspond to HERMES data [12] and the 
open points are obtained by HARDPING 2.0 
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